Minimal relativity and 3 Si- 3 D\ pairing in symmetric nuclear matter 
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We present solutions of the coupled, non-relativistic 3 Si- 3 Dj gap equations for neutron-proton 
pairing in symmetric nuclear matter, and estimate relativistic effects by solving the same gap 
equations modified according to minimal relativity and using single-particle energies from a Dirac- 
Brueckner-Hartree-Fock calculation. As a main result we find that relativistic effects decrease the 
value of the gap at the saturation density kp = f .36 fm _1 considerably, in conformity with the lack 
of evidence for strong neutron-proton pairing in finite nuclei. 
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The size of the neutron-proton (np) 3 Si- 3 Di energy 
gap in symmetric nuclear matter has been a much de- 
bated issue since the first calculations of this quantity ap- 
peared. While solutions of the BCS equations with bare 
nucleon-nucleon (NN) forces give a large energy gap of 
several MeV's at the saturation density kp = 1-36 fm _1 
(p = 0.17 fin -3 ) there is little empirical evidence 

from finite nuclei for such strong np pairing correlations. 
One possible resolution of this problem lies in the fact 
that all these calculations have neglected contributions 
from the so-called induced interaction. Fluctuations in 
the isospin and the spin-isospin channel will probably 
make the pairing interaction more repulsive, leading to a 
substantially lower energy gap ||. One often neglected 
aspect is that all non-relativistic calculations of the nu- 
clear matter equation of state (EOS) with two-body NN 
forces fitted to scattering data fail to reproduce the em- 
pirical saturation point, seemingly regardless of the so- 
phistication of the many-body scheme employed. For ex- 
ample, a Brueckner-Hartree-Fock (BHF) calculation of 
the EOS with one of the Bonn potentials would typi- 
cally give saturation at hp = 1.6-1.8 fm — . In a non- 
relativistic approach it seems necessary to invoke three- 
body forces to obtain saturation at the empirical equilib- 
rium density. This leads one to be cautious when talk- 
ing about pairing at the empirical nuclear matter satura- 
tion density when the energy gap is calculated within a 
pure two-body force model, as this density will be below 
the calculated saturation density for this two-body force, 
and thus one is calculating the gap at a density where 
the system is theoretically unstable. One even runs the 
risk, as pointed out in Ref. [bj, that the compressibil- 
ity is negative at the empirical saturation density, which 
means that the system is unstable against collapse into 
a non-homogeneous phase. A three-body force need not 
have dramatic consequences for pairing, which after all 
is a two-body phenomenon, but still it would be of in- 
terest to know what the 3 iSi- 3 .Di gap is in a model in 
which the saturation properties of nuclear matter are re- 
produced. If one abandons a non-relativistic description, 
the empirical saturation point can be obtained within 
the Dirac-Brueckner-Hartree-Fock (DBHF) approach, as 
first pointed out by Brockmann and Machleidt (7). This 



might be fortuitous, since, among other things, impor- 
tant many-body effects are neglected in the DBHF ap- 
proach. Nevertheless, we found it interesting to inves- 
tigate 3 S\- 3 Di pairing in this model and compare our 
results with a corresponding non-relativistic calculation. 

The first ingredient in our calculation is the self- 
consistent evaluation of single-particle energies in sym- 
metric nuclear matter starting from the meson-exchange 
potential models of Machleidt and co-workers J8|. For 
the non-relativistic (NR) calculations we use the BHF 
method, while the DBHF scheme is used in the relativis- 
tic (R) calculation. Details of both approaches are found 
in Refs. Since BHF and DBHF are computa- 

tionally very similar, we will here content ourselves with 
giving a brief description of the DBHF method. In this 
scheme, the single-particle energies and binding energy 
of nuclear matter is obtained using a medium renormal- 
ized NN potential G denned through the solution of the 
G-matrix equation 



G(u) = V + VQ- 
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where u> is the unperturbed energy of the interacting nu- 
cleons, V is the free NN potential, Hq is the unperturbed 
energy of the intermediate scattering states, and Q is 
the Pauli operator preventing scattering into occupied 
states. Only ladder diagrams with two-particle interme- 
diate states are included in Eq. ([!]). In this work we solve 
Eq. (jl]) using the Bonn A potential defined in Table A. 2 
of Re f. M . This potential model employs the Thomp- 
son ||7|jll[ reduction of the Bethe-Salpeter equation, and 
is tailored for relativistic nuclear structure calculations. 
For the non-relativistic calculation we employ the Bonn 
A potential with parameters from Table A.l in Ref. [pj. 
This model employs the Blankenbecler-Sugar (BbS) re- 
duction of the Bethe-Salpeter equation, and is therefore 
suited for non-relativistic calculations. For further de- 
tails, see Refs. @,|,@. 

The DBHF is a variational procedure where the single- 
particle energies are obtained through an iterative self- 
consistency scheme. To obtain the relativistic single- 
particle energies we solve the Dirac equation for a nucleon 
in the nuclear medium, with c = h = 1, 
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m + Y,[p))u{p, s) = 0, 



(2) 



where m is the free nucleon mass and u(p, s) is the Dirac 
spinor for positive energy solutions, p = (p , p) being 
a four momentum and s the spin projection. The self- 
energy S(p) for nucleons can be written as 

£(p) = £ s (p)- 7o I] ( V)+7-p£ y (p). (3) 

Since S v << 1 |?||l^], we approximate the self-energy by 

£ « E s - 7o E° = C/ s + CV, (4) 

where L^s is an attractive scalar field and Uv is the time- 
like component of a repulsive vector field. The Dirac 
spinor then reads 



u(p, s) 



' E p + rh 
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(5) 



Ep+ih 



where Xs is the Pauli spinor and terms with tilde like 
E p = y/p 2 + fh? represent medium modified quantities. 
Here we have defined [f],[l2| m = m + Us- The single- 
particle energies s p can then be written as 



(p| 7 ' P + m \p) + ' 



E„ 



U v , 



(6) 



where the single-particle potential w p is given by u p — 
Usm/Ep + Uy and can in turn be defined in terms of the 
G-matrix 



p'<k F P • 
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where p, p' represent quantum numbers like momentum, 
spin, isospin projection etc of the different single-particle 
states and kp is the Fermi momentum. Eqs. are 
solved self-consistently starting with adequate values for 
the scalar and vector components Us and Uv- The en- 
ergy per particle can then be calculated from 
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where the subscripts and 2 denote S and -D states re- 
spectively, Vu> is the free momentum-space NN interac- 
tion in the relevant channel, Ao and A 2 are the S and D 
state gap functions respectively, and E(k) is the quasi- 
particle energy given by 



E(k) = V(£fe-/i) 2 + A (fc)2 + A 2 (fc)2, (11) 
where fi is the chemical potential. The quantity 



D F = v/A (fc F ) 2 + A 2 (fc F ) 2 



(12) 



will in the following be referred to as the energy gap, 
according to the conventional definition [Q-Q. For 3 Si- 
3 D\ pairing it is also necessary to solve the equation for 
particle number conservation 
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In Fig. ^ we show the EOS obtained in our non- 
relativistic and relativistic calculations. The non- 
relativistic one fails to meet the empirical data, while 
the relativistic calculation very nearly succeeds. 

Having obtained in-medium single-particle energies, 
we proceed to solve the coupled gap equations for 3 Si- 
3 Di pairing. Employing an angle-average approximation, 
these can be written || 



(9) 



for /i self-consistently together with Eqs. (||) and (|To|). 

In the non-relativistic calculation, we have used the 
Bonn A potential with parameters from table A.l in Ref. 
||. The results are shown in Fig. |[ We found a large 
energy gap at the empirical saturation density, around 
6 MeV at kp — 1.36 fm _1 , in agreement with earlier 
non-relativistic calculations 

Recently, several groups have developed relativistic for- 
mulations of pairing in nuclear matter |l3|-[l5f| , and have 
applied them to Sq pairing. The models are of the 
Walecka-type |lj] in the sense that meson masses and 
coupling constants are fitted so that the mean-field EOS 
of nuclear matter meets the empirical data. In this way, 
however, the relation of the models to free-space NN scat- 
tering becomes somewhat unclear. An interesting result 
found in Refs. |l^-p^|| is that the 1 So energy gap vanishes 
at densities slightly below the empirical saturation den- 
sity. This is in contrast with non-relativistic calculations 
which generally give a relatively small, but non- vanishing 
1 5o gap at this density, see for instance |l6|-[l9|. 

With these results in mind, we found it interesting to 
consider relativistic effects on 3 5'i- 3 Di pairing in nuclear 
matter, which, to our knowledge, has not been done be- 
fore. A simple way of doing this is to incorporate minimal 
relativity in the gap equation, thus using DBHF single- 
particle energies in the energy denominators and modi- 
fying the free NN interaction by a factor m 2 /EkEk> [p0| . 
With this prescription, we obtained the results shown in 
Fig. H (short-dashed line). As can be seen, the gap at 
the empirical saturation density is reduced from 6 MeV 
to nearly zero. 
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Let us try to understand the difference between the rel- 
ativistic and the non-relativistic calculation. First of all, 
we point out that it is well known that the introduction 
of relativity in the many-body problem leads to increased 
repulsion at densities at and above the saturation den- 
sity Even a slight increase in the repulsion might 
have large consequences for the energy gap, as the gap 
depends exponentially on the interaction at the Fermi 
surface. One can obtain a numerical estimate of the ef- 
fect as follows. If one takes the weak-coupling limit of 
the gap equations, Eqs. (^) and (|Io|), it is easy to show 
that one obtains the same form for the energy gap as for 
1 S pairing, 



Dp = 26 e exp 



N(k F )V pal 



(14) 



where Se is an appropriate energy interval, N(k F ) = 
rim* kp /2tt 2 % 2 is the density of states at the Fermi sur- 
face, to* is the nucleon effective mass, and f2 is the volume 
occupied by the system. The pairing interaction at the 
Fermi surface, V pa i r , is given by 



V, 



y/(Vss~V DD y+W s 2 D - V ss ~ V SD 



pair 



, (15) 



where Vss = V m {k F ,k F ), V S d = V m {k F ,k F ) = 
V2o(kp,k F ), and Vdd = V22(k F ,k F ). With our rela- 
tivistic approach to the gap equation, the correspond- 
ing weak-coupling expression for the gap is obtained by 
replacing V pa i r with rr? j E\ F V vair and using the rela- 
tivistic expression for the density of states instead of 
the non-relativistic one. We will now consider the sat- 
uration density, and take k F « 1.4 £m _ . The non- 
relativistic single-particle spectrum was parameterized 
as €k = k 2 /2m* + U . At k F = 1.4 fnV 1 the values 
of the parameters were m* jm — 0.6751, Uq — —97.2755 
MeV. For the relativistic single-particle spectrum the rel- 
evant quantities were U s = -384.89 MeV, U v = 300.18 
MeV. We found N R (k F = 1.4 rmfm^ 1 ) / N NR {k F = 
1.4 fm" 1 ) ps 1, and then 



exp 



1 



N NR {lA)V pa 

1/4 



D NR 

2Se 



where the superscript R refers to relativistic quantities, 
NR to non-relativistic ones. If we make the common 
choice Se = €p R and use Dp m 6 MeV, we obtain 



D NR 
F 



0.5, 



thus, the introduction of relativity in the gap equation 
suppresses the gap at the saturation density by a fac- 
tor of two. This argument makes it reasonable that rel- 
ativistic effects reduce the gap. That the reduction is 



larger in the full calculation than in this simple estimate 
is understandable, since we in the weak coupling approx- 
imation neglect the momentum-dependence of the inter- 
action. More specifically, the repulsive high-momentum 
components are left out, and these will reduce the gap 
further. 

It is also interesting to obtain the ratio D R /Dp R at 
the respective saturation densities for the relativistic and 
non-relativistic calculations. The non-relativistic EOS 
saturates at k F 1.8 fm" 1 . At this density, we had nu- 
merical problems with solving the gap equations, some- 
thing which may occur when the gap is small. However, 
using the weak-coupling expression for the gap, we could 
estimate D§(kp = 1.4 fm^ 1 )/D^ R (k F = 1.8 fnV 1 ). 
First we used the non-relativistic gaps in the density 
range k F = 1.2-1.4 fm -1 to calculate NN R (k F )V pa i r (k F ), 
and fitted the results with a quadratic polynomial in kp. 
From this fit we estimated N^ R (k F = l.8)V pa i r (kp — 
1.8) « 0.257. Then 



N NR {l.8)V patr (l.8) _ 0.257 



N NR (lA)V pair (lA) 
=S> V pair (l.8) 



0.333 

iVjvfl(1.4) 0.257 
N NR (1.8) 0.333 
0m5V pair {lA) 



V p air (1-4) 



where we have used m*(1.4)/m = 0.675, m*(1.8)/m = 
0.5834 in the densities of states. We then formed the 
ratio D R (1A)/D F ' R {1.8) and obtained 



n R 

U F 
D NR 
F 



= exp[ 



1 



N NR (l.8)V pair (lA) 
X NR (1.8) £|(1.4) 1 

in 2 



N R (1A) m 2 (lA) 0.695 



since we found N NR (1.8)/N R (1A) « 1.15, 
E% (1.4) /to 2 (1.4) w 1.25, thus making the expression 
in the inner parenthesis « 0. Although this argument is 
only indicative, it makes it reasonable to assume that the 
non-relativistic gap will be very small at the calculated 
non-relativistic saturation density. 

In this Letter, we have presented non-relativistic and 
relativistic calculations of the 3 S'i- 3 Di np gap in sym- 
metric nuclear matter. The non-relativistic calculations 
gives a large gap of approximately 6 MeV at the empirical 
saturation density. In the relativistic calculation we find 
that the gap is vanishingly small at this density. This 
is the main result of this Letter. Non-relativistic calcu- 
lations with two-body interactions will in general give a 
saturation density which is too high, an example of which 
is shown in Fig. |l|. Thus in the non-relativistic approach 
we are actually calculating the gap at density below the 
theoretical saturation density, and one may question the 
physical relevance of a large gap at a density where the 
system is theoretically unstable. If one looks at the gap 
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at the calculated saturation density, it is in fact close to 
zero. In the DBHF calculation we come very close to 
reproducing the empirical saturation density and bind- 
ing energy, and when this is used as a starting point for 
a BCS calculation, we find that the gap vanishes, both 
at the empirical and the calculated saturation density. 
That the DBHF calculation meets the empirical points 
is perhaps fortuitous, as important many-body diagrams 
are neglected and only medium modifications of the nu- 
cleon mass are accounted for. However, a recent inves- 
tigation where medium modifications of meson masses 
were included, showed that the results of the DBHF do 
not change very much, and in particular the saturation 
properties are still very good p^] . Nevertheless, the es- 
sential property which is needed in all non-rclativistic 
models to get to the empirical point is an increased re- 
pulsion at and around the empirical saturation density. 
Regardless of the mechanism, this may reduce the pair- 
ing gap dramatically. The main point we wish to make in 
this Letter is thus that the inclusion of these additional 
repulsive effects may destroy pairing at the empirical sat- 
uration density, a result which is in accordance with the 
lack of empirical evidence for the strong np pairing pro- 
duced in non-relativistic calculations. 
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FIG. 1. EOS for symmetric nuclear matter with the NN potentials and many-body methods described in the text. 
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FIG. 2. 3 Si- 3 D! energy gap in nuclear matter, calculated in non-relativistic (full and long-dashed lines) and relativistic 
(short-dashed line) approaches. 
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